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The possibility to use 7-ray data from the Galactic Center (GC) to constrain the cosmological 
evolution of the Universe in a phase prior to primordial nucleosyntesis, namely around the time of 
cold dark matter (CDM) decoupling, is analyzed. The basic idea is that in a modified cosmological 
scenario, where the Hubble expansion rate is enhanced with respect to the standard case, the CDM 
decoupling is anticipated and the relic abundance of a given dark matter (DM) candidate enhanced. 
This implies that the present amount of CDM in the Universe may be explained by a Weakly 
Interacting Massive Particle (WIMP) which possesses annihilation cross section (much) larger than 
in standard cosmology. This enhanced annihilation implies larger fluxes of indirect detection signals 
of CDM. We show that the HESS measurements can set bounds for WIMPs heavier than a few 
hundreds of GeV, depending on the actual DM halo profile. These results are complementary to 
those obtained in a previous analysis based on cosmic antiprotons. For a Moore DM profile, 7-ray 
data limit the maximal Hubble rate enhancement to be below a factor of 100. Moreover, a WIMP 
heavier than 1 TeV is not compatible with a cosmological scenario with enhanced expansion rate 
prior to Big Bang Nucleosynthesis (BBN). Less steep DM profiles provide less stringent bounds, 
depending of the cosmological scenario. 

PACS numbers: 95.35.+d,95.36.+x,98.80.-k,04.50.+h,96.50.S-,98.70.Sa,98.80.Cq 



I. INTRODUCTION 



In a recent paper |l[ we discussed the possibility to de- 
rive limits to the expansion rate of the Universe around 
the time when cold dark matter (CDM) decouples from 
the thermal bath, by using the fact that today these dark 
matter (DM) particles form the halo of our Galaxy and 
they may annihilate producing antiprotons. The mea- 
sured antiproton flux has been shown to be compatible 
with the expected background originated by standard 
cosmic-ray spallation: this fact leads to the possibility 
to use antiprotons as a powerful tool for constraining 
DM properties (see e.g. Refs. By assuming that 

DM is in the form of a generic Weakly Interacting Mas- 
sive Particle (WIMP) candidate, we obtained relatively 
stringent bound on the Universe dynamics in a period 



prior to the primordial nucleosynthesis phase (preBBN), 
which is not directly constrained by other observations. 
We showed that, despite the large uncertainties in the 
knowledge of the galactic propagation of antiprotons 0] , 
bounds on the Hubble rate enhancement ranging from a 
factor of a few to a factor of 30 are present for DM masses 
lighter than 100 GcV, while for a mass of 500 GeV the 
bound falls in the range 50-500. These bounds loosen for 
heavier DM particles. 

We remind here that the possibility to set a bound 
to the enhancement of the Hubble expansion rate in the 
early Universe by means of WIMP indirect detection sig- 
nals relies on the fact that a larger Hubble rate induces 
an anticipated DM decoupling and an ensuing larger relic 
abundance, for a fixed annihilation cross section. This 
may happen, for instance, in scalar-tensor cosmology 0], 



2 



in quintessence models with a kination phase 0] , and for 
anisotropic expansion and other models of modified ex- 
pansion @, 0). 

When the decoupling is anticipated, a DM relic abun- 
dance able to explain the current observational data is 
obtained for larger DM annihilation cross sections, as 
compared to the standard cosmological case. This means 
that larger indirect detection signals (which are propor- 
tional to the DM annihilation cross section) are pre- 
dicted: the larger the enhancement of the Hubble rate, 
the larger the indirect detection signal [if. Comparison 
with experimental data reflects in limits on the DM pos- 
sible Hubble rate enhancement. The actual results de- 
pend also on the specific cosmological model: in Ref. 
|l| we considered a general parametrization of the Hub- 
ble rate temperature dependence, and for dcfinitcness we 
studied different classes of models: a Randall-Sundrum 
branc cosmology scenario 0], a kination scenario 0], a 
scalar-tensor cosmology scenario and a simple case 
where the Hubble rate is just boosted by a constant fac- 
tor. The same modcllizations will be considered in the 
present study. 

In this paper we extend our analysis by including also 
the possibility to use the data on 7-rays coming from the 
GC. DM particles may annihilate also to photons and 
the current experimental observations on the 7-ray flux, 
especially the ones coming from the GC where a large 
concentration of DM is expected, may be useful for our 
scope. With 7-rays we have to face the large uncertainty 
coming from the DM halo profile toward the GC, which 
is largely unknown. Nevertheless we will show that for 
typical DM halo profiles, bounds to the preBBN Hubble 
rate may be set. 

The experimental data we will use are those from the 
GC obtained by the HESS and EGRET [ll| tele- 

scopes. Since the HESS data refer to large energies 
(above few hundreds of GeV), we will be able to derive 
bounds for large-mass DM (above few hundreds of GeV): 
this nicely complements our previous results on antipro- 
tons which were quite stringent for low mass WIMPs, 
but loose for heavy ones. We will also show that EGRET 
data in the GeV region are able to set limits in the inter- 
mediate WIMP mass range (at least for the most steep 
DM profiles). 

Our analysis will be independent of the specific DM 
candidate: the only assumption is that the WIMP we are 
considering is responsible for the observed DM amount 



0.092 < n CDM /i 2 < 0.124 



(1) 



where ficDM is the density parameter of CDM, and its 
annihilation cross section is dominantly temperature- 
independent (or s-wave). Deviations from this situation 
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FIG. 1: Different models for the Hubble rate enhancement 
function A(T), defined in Eq. @. The Figure is taken from 
Ref. Q. Notice that the evolution of the Universe runs from 
right to left. The solid (red) line has a slope parameter v = 
— 1, the dashed (green) line v — 0, the dash-dotted (blue) 
v = 1 and the dotted (purple) line has v = 2. All curves refer 
to 77 = 10 2 , T rc = 10~ 2 GeV and T t = 17.3 GeV (this freeze- 
out temperature refers, e.g., to a particle with mass of 500 
GeV and annihilation cross section (<r ann u ) = Kr 7 GeV~ 2 ). 



will lead to changes similar to those discussed in Ref. [l| 
for the antiproton case, to which we refer for additional 
discussion. 



II. MODELS WITH INCREASED PRE-BBN 
EXPANSION 

As mentioned in Sect. HI it is a common feature of some 
cosmological models to predict that the expansion rate 
H(T) in the early Universe is larger than the Hubble 
expansion rate Hqr(T) in standard cosmology. Quite 
generically, we may introduce a function A(T) to quantify 
the enhancement of the Hubble rate: 

H(T) = A(T)H GR (T) at early times (2) 

H(T) = H GK (T) at later times. (3) 

In Ref. [l[ we introduced a parametrization of the en- 
hancement function A{T) which was shown to be appli- 
cable for important models like some scalar-tensor grav- 
ity models, some models with a kination phase and also 
some specific brane-world model. In this paper, we are 
going to consider the same parametrization, i.e.: 



(| 



tanh 



T - T r< 



(4) 
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for temperatures T > T IC and A(T) = 1 for T < T rc . 
The enhancement function is shown in Fig. [1] for some 
specific parameter choices. The hyperbolic tangent serves 
to assure that A(T) goes continuously to "1" , and H — > 
i?GR 5 before some "re- entering" temperature. T ro . We 
must require T rc > 1 MeV to make sure not to be in 
conflict with the predictions of BBN and the formation of 
the Cosmic Microwave Background (CMB). For T > T IC 
(and rj ^> 1) we have approximately that: 



(5) 



Thus, Tf is the normalization temperature at which 
A(Tf) = r/. As in Ref. Q we take Tf to be the tem- 
perature at which the WIMP DM candidate freezes out 
in standard cosmology. This means that rj, as defined in 
our paramctrization, is the enhancement of the Hubble 
rate at the time of the WIMP freeze-out. We will de- 
rive our results as bounds on rj for different cosmological 
models, characterized by the temperature-evolutionary 
parameter v. v = 2 refers to the Hubble rate evolution 
in a Randall-Sundrum type II brane cosmology scenario 
of Ref. 0; v = 1 is the typical kination evolution, dis- 
cussed e.g. in Ref. Q); f = — 1 is representative of the 
behavior found in scalar-tensor cosmology in Ref. 0]. 
The trivial case v = refers to an overall boost of the 
Hubble rate. 

For more details on the modified cosmological scenar- 
ios, the calculation of the relic abundance in these mod- 
els, including some analytical results and discussion, we 
refer the reader to Ref. [l|. 



III. THE 7-RAY SIGNAL FROM THE GC 

The most recent observations of the GC have been per- 
formed by the HESS Collaboration in the hundreds of 
GeV - few tens of TeV energy range. In Refs. 0, E3, 03] 
they have reported on the spectrum of very high energy 
7-rays from a point-like source in the GC, with an un- 
precedented spatial resolution, going down to a solid an- 
gle of about 10~ 5 sr. The HESS spectrum of the central 
7-ray source exhibits a clear power-law shape, with a 
spectral index of 2.20-2.25[I3|- Diffuse 7-ray emission 
extended along the galactic plane has been reported by 
the same Collaboration in Ref. 



14[ . The measured spec- 



trum follows a power-law with spectral index near to 2.30 
and has been shown to be compatible with a source of 
locally accelerated protons interacting with giant molec- 
ular clouds which are extended both in longitude and in 
latitude. This diffused component contributes to the cen- 
tral source only for a small fraction (10-15%) 14 1. 
The GC hosts more than one potential sources of 7-rays, 



whose nature is still not clear. The most motivated as- 
trophysical sources rely on particle acceleration near the 
supermassive black hole Sgr A* located at the center of 
our Galaxy, or in the region of the supernova remnant Sgr 
A. The 7 radiation is produced from accelerated charged 
particles (mostly protons) interacting with the ambient 
matter or radiation. Another intriguing possibility re- 
sides in the 7-ray emission resulting from the annihilation 
of DM particles. Cosmological simulations of hierarchical 
structure formation predict a significant density cusp in 
the central parts of the galaxies. In that region annihi- 
lation of DM particles would be strongly enhanced, with 
extraordinary high expected fluxes for the annihilation 
products, such as 7-rays. The angular region explored 
by HESS is of the same order of the one of the probable 
black hole, or of the DM cusp. 

We consider here a generic WIMP which composes the 
galactic DM. The flux <I> 7 (i? 7 , ip) of 7-rays of energy E 1 
originated from the WIMP pair annihilation and coming 
from the angular direction ip is given by: 



1 (av) dN y 1 , , N 



(6) 



where {crv}o is the present annihilation cross section times 
the relative velocity averaged over the galactic velocity 
distribution function. dN^/dEy is the energy spectrum 
of 7-rays originated from a single WIMP pair annihila- 
tion and has been calculated by means of a Monte Carlo 
simulations with the PYTHIA package [3] as described 
in Ref. [ljij]. For definiteness, as we have done also for 
the antiproton analysis of Ref. [l| , we are assuming the 
7-ray energy spectrum originated by a bb quark pair. 
A different annihilation final state will not change sub- 
stantially our results, much less than the astrophysical 
uncertainties. 

The quantity I(ip) is the contribution of the squared 
DM density distribution along the line of sight (l.o.s.): 



(7) 



l.o.s 



Here i/j is the angle between the l.o.s. and the line point- 
ing toward the GC (cos ip = cos I cos b, I and b being the 
galactic longitude and latitude, respectively). If Eq.j7|) 
is used for comparison with experimental data, it must 
be averaged over the telescope observing angle A?/;: 



(8) 



IatP = -r— - / dip. 



The geometric factor I(ip) depends quadratically on 
the DM density profile and is very sensitive to its fea- 
tures, especially in the GC region where predictions 
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Isothermal 


NFW 


log-slope 


a = 1.2 


Moore 


18.9 


6892 


10229 


98743 


7.7 • 10 6 



TABLE I: Values for Javj in Eq.® (in units of GeV 2 cm 6 kpc). See text for details. 



for p(r) differ mostly. The most common spherically- 
averaged density profiles can be parametrized as: 



P( r ) =Pl[ — 
r 



7 r 



1 + (R Q /a)° 
1 + {r/a) a 



(9) 



where r = \f\, Rq = 8 kpc is the distance of the So- 
lar System from the GC along the galactic plane, a is a 
scale length and pi is the total local (solar neighborhood) 
DM density. In particular, the isothermal, cored density 
profile is obtained with (a,f3, 7) = (2.2,0), the Navarro, 
Frcnk and White (NFW) profile with (a, (3, 7) = 
(1,3,1) and the Moore ct al. profile [lj| with (a,/3, 7) 
= (1.5,3,1.5). We also consider a non-singular DM den- 
sity distribution function derived from extensive numer- 
ical simulations, whose asymptotic regime is well fitted 
by a logarithmic slope [19( : 



p(r) = P-2 cxp <; 

a 



r 
r-2 



1 



(10) 



where r_2 is the radius where the logarithmic slope is 
5 = —2, and p-2 = p(f-2)- The DM density predicted 
by these profiles at the GC (for very small r) varies so 
strongly that also the predicted DM signals may differ by 
several orders of magnitude. The profile as steep as 1.5 
is disfavored by the most recent numerical simulations, 
which seem to indicate a power law index not exceeding 
1.2 20|, |2l|. We notice, however, that the experiments 



considered in this paper have a spatial resolution which 
is much narrower that the typical resolution size of nu- 
merical simulations. The simulated DM densities in the 
GCs are thus mere extrapolations. Moreover, we must be 
aware of the fact that the results from many rotational 
curves for galaxies of different morphological types are 
hardly explained by central DM cusps. Instead, they are 
more easily fitted by cored DM distributions, flattened 
towards the central region of the galaxy. 

In Table [J we show the values of the geometrical factor 
Ia,p for the HESS telescope aperture (Aip = 10~ 5 sr) 
and for various DM density profiles. The first column 
refers to an isothermal density distribution with a core 
a = 3.5 kpc and the second to a NFW with a scale length 
a = 25 kpc. The third is for the log-slope of Eq. (fTOl) 
with the parameters of the distribution Gl in Ref. [l 91 ] : 
a = 0.142, r_ 2 = 26.4 kpc and p 2 = 0.035 GeV cm" 3 . 
The fourth column is the result for a profile obtained 
with (a,/3,7) = (1.2,3,1) with a=25 kpc 0,110] and the 
last column refers to a Moore et al. profile with scale 



length a = 30 kpc. The value of pi can be determined for 
each density profile requiring the compatibility with the 
measurements of rotational curves and the total mass of 
the Galaxy 22j. For definiteness, we have fixed pi = 0.3 
GeV cm~ 3 for all the density profiles in Table HI We 
notice that the parameter pi enters as a mere scaling 
factor in Eq. [6j the effect of varying pi is therefore easily 
taken into account. 



IV. CONSTRAINING THE HUBBLE RATE 
WITH 7-RAYS 

In order to derive the constraint on the enhancement of 
the Hubble rate, we first find the upper bound on {av)o 
determined by the 7-ray observations, for any given halo 
profile. 

Let us first explain our analysis for the HESS obser- 
vation of high energy 7-rays from the GC. Here we have 
added the expected WIMP signal to a background that 
follows a power-law, kE~ r . The normalization and the 
index of the power-law are taken as free parameters. For 
each point in the grid-scan of (k, T) we find the maxi- 
mum allowed value of {cfv)q, which we statistically define 
as the value of the cross section at which the reduced \ 2 
equals 3. The reduced \ 2 is defined as: 



A, red 



AT ^ 



N 



$ x (n) + ^ 5 (n)-$o bs (n) 
cobs {n) 



(11) 



where N is the number of degrees of freedom, & x (n) is the 
expected WIMP flux, as calculated using Eq. ©, at the 
energy of the n-th data-point ($ x (n) obviously vanishes 
when E n > m x ). The background flux is assumed to be 
$s(n) = kE~ r , while the observed flux is <E>obs(«-)- The 
1-ct error of the observation is denoted uobs (n) . For the 
analysis, we use the 17 datapoints from the HESS 2003 
observation of the GC [^, [lj]. They cover the energy 
range 280 GeV - 8.83 TeV. These data can therefore be 
used to constrain WIMP heavier than about 300 GeV 
(the annihilation process in the Galaxy occurs almost at 
rest: therefore there is a kinematic cut off at the WIMP 
mass for the 7-ray energy). 

The x 2 analysis, as described above, gives us the upper 
limit of the cross section for each point in the (k, T) grid 
for fixed WIMP mass and halo profile. Finally we extract 
the grid-point which gives the biggest value of the cross 
section. This upper bound on (crv)o is shown in Fig. [5J as 
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FIG. 2: Bounds on the WIMP annihilation cross section (ov)q 
as a function of the WIMP mass. The curved lines at large 
WIMP mass show the upper bounds derived from the HESS 
observation of 7-rays from the GC. The derivation has been 
made for five different DM halo profiles. From top to bottom 
these are: the isothermal model, the NFW model, the 'log' 
slope, a power-law slope with index 1.2, the Moore et al. 
profile (see text for more details). The dashed lines show 
analogous upper limits derived using the 7-ray data from the 
EGRET detector, for the same set of galactic halo models. 
The EGRET limits are plotted only in the mass interval which 
is relevant for the analysis of this paper. The slanted solid 
lines show the upper limits coming from the observations of 
cosmic antiprotons Q|. The central line refers to the best 
estimate for the antiproton DM signal. The upper and lower 
lines refer to the astrophysical uncertainties in the galactic 
propagation parameters [l|, 0|. Finally, the horizontal dot- 
dashed lines shows the lower bound on the cross section as 
derived from the WIMP relic density constraint assuming the 
standard cosmological model and a temperature independent 
WIMP annihilation cross section. 



a function of the WIMP mass for the five different halo 
profiles discussed above. As expected, the result depends 
strongly on the halo profile. 

Before we continue to explain how the upper limit on 
the cross section was derived from other observational 
data, let us show some examples of the differential pho- 
ton production which correspond to the upper limit on 
the WIMP annihilation cross section. In Fig. [3] we 
show the result for the NFW halo profile and for two 
different masses. The flux is calculated using Eq. (J6]) 
and inserting the upper bound on (eru)o for the given 
mass and halo profile. Also shown are the HESS 2003 
observations of the GC as well as the fitted power-law 
background. Note that as the parameters of the back- 
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FIG. 3: 7-ray spectra as a function of energy. The data 
points refer to HESS 0, [l3|. The solid and dashed (green) 
curves refer to the maximal allowable contribution to the 7- 
ray flux from WIMP annihilation: the cases of 500 GeV and 5 
TeV WIMPs are plotted. The solid and dashed (red) straight 
lines refer to a standard-source power-law contribution to the 
HESS data, as obtained by our fit (see text for more details). 



ground are treated as free parameters, they are differ- 
ent for the different WIMP masses and halo profiles. In 
the shown examples, the background parameters associ- 
ated with the upper bound on the WIMP cross section 
are (fc,T) = (0.177 ■ 10" 6 cm" 2 s" 1 TeV" 1 sr~ 1 ,2.00) 
and (k,T) = (0.133 • 10" 6 cm" 2 s" 1 TeV" 1 sr" 1 , 1.83) 
for WIMP masses of 500 GeV and 5 TeV respectively. 



Coming back to the limit on the WIMP annihilation 
cross section, Fig. [2]shows also the upper limit as derived 
from the EGRET. The EGRET data [ill s P an from en- 
ergies of around 0.039 GeV to around 14.9 GeV with an 
angular resolution given by the longitude-latitude aper- 
ture: |A/| < 5°, |A6| < 2°. The geometric factor for 
the EGRET experiment has been taken from Ref. [l6l ]. 
for the same density distribution functions described in 
the previous Section. At these energies there is a 7- 
ray background from nucleonic reactions between cosmic 
rays and interstellar medium, electron brcmsstrahlung 
and inverse Compton. We assume it to be the same as 
the one described in Ref. We apply the following 

analysis only for the two highest EGRET energy bins, 
i.e. E n ~ 6.2 GeV and E n ~ 14.9 GeV, which are the 
most constraining for the masses we are dealing with. 
For a given WIMP mass and halo profile we make a scan 
in the WIMP annihilation cross section to find its upper 
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FIG. 4: Summary of the DM indirect detection limits on the 
WIMP annihilation cross section (for a temperature indepen- 
dent cross section). The allowed region lies above the hori- 
zontal dot-dashed lines (which refer to the relic density con- 
straint) and below the slanted/curved lines, for any given DM 
density profile. The slanted part of the upper bound is due 
to cosmic antiprotons, while the curved part to 7-rays from 
the GC (from HESS data at large WIMP masses, EGRET at 
intermediate masses). 

limit taken to be the 2tr bound: 

$ x (n) + $s(n)-$ bsW <2a ohs (n) (12) 

The observational error, that we use, only includes the 
statistical error. We use the limit from the most con- 
straining of the two data points. For the mass range rel- 
evant here, it is always the highest energy bin which pro- 
vides the bound, except for a WIMP mass of 100 GeV. In 
Fig. [2] we show only the WIMP range where the EGRET 
data could be of importance for the analysis of this paper. 
Finally, Fig. [2] shows, as slanted solid lines, the upper 
limit on (ctv)q derived in our previous analysis 0] of the 
observational data on cosmic antiprotons. The central 
line refers to the best estimate for the antiproton DM 
signal, i.e when the best set of astrophysical parameters 
are used in the calculation of the diffusion processes for 
the galactic cosmic rays 0]- The upper and lower lines 
refer to the uncertainty band arising from astrophysi- 
cal uncertainties in the galactic propagation parameters 
[H, 0] ■ The horizontal dot-dashed line denote instead the 
lower limit on the annihilation cross section derived from 
the cosmological bound on the WIMP relic density Eq. 
([1]) in standard cosmology. 

The analysis of Fig. [2] shows interesting properties. In 
addition to the already discussed bound on (o~v)q from 



antiprotons, which sizeably constrains the maximal al- 
lowable annihilation cross section, especially for WIMPs 
lighter than a few hundreds of GeV, we now also have the 
bounds coming from 7-rays from the GC, which instead 
are relevant for heavy WIMPs. This is a consequence of 
the fact that HESS data refer to an energy range from a 
few hundreds of GeV up to few TeV. The figure clearly 
shows that, in order to have bounds on (o~v)o, the sig- 
nal must be quite sizeable and this is possible only for 
very steep DM profiles like the Moore and NFW ones. 
In the case of the Moore profile, we have a tension be- 
tween the 7-rays HESS data and the cosmological limit 
even for standard cosmology (a situation analogous to 
the one already observed for the antiproton signal pro- 
duced by light WIMPs UM)- Should the DM profile 
be the Moore one, very little room would be allowed for 
(crw)o- This bound would also imply a finite possible 
range for the WIMP mass: from 50-100 GeV to about 1 
TeV (even for standard cosmology), as a combination of 
cosmological data on the amount of DM, the antiproton 
component of cosmic-rays and the 7-ray signal for the 
GC. Less steep profiles are clearly much less constrain- 
ing, and the isothermal case is always ineffective, being 
always less relevant than the antiproton bound. We re- 
mind that since antiprotons diffuse in the Galaxy, their 
flux is only mildly dependent on the DM profile [2j]: the 
antiproton bounds arc therefore practically unaffected by 
the choice of different halo shapes. 

A summary of the different bounds derived in our anal- 
ysis is given in Fig. [?J where we combine the upper limits 
on the WIMP annihilation cross section from the obser- 
vations of the cosmic antiprotons and of the GC 7-rays, 
as observed by both EGRET and HESS. For any mass 
and halo profile we take the bound which is most con- 
straining. The combined upper limit can be seen in Fig. 
2] for three halo profiles. 



V. THE MAXIMAL ENHANCEMENT OF THE 
HUBBLE RATE 

An increased Hubble rate in the early Universe would 
increase the relic density of a given WIMP as compared 
to the situation in standard cosmology. The increase of 
the relic density is due to an anticipated freeze-out of the 
WIMP, as the annihilation rate cannot keep up with the 
expansion rate as long as in the standard case. WIMPs 
which satisfy the density constraint of Eq. (JT|) in the mod- 
ified scenario would thus be undcrabundant in standard 
cosmology. Because of the inverse proportionality be- 
tween the WIMP relic density and annihilation cross sec- 
tion, the WIMPs which fulfill the density constraint in 
the modified cosmologies have larger cross section than 
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FIG. 5: Upper bound on the Hubble rate enhancement- 
parameter r\. The bound is shown as a function of the WIMP 
mass and has been derived by combining the constraints on 
the WIMP relic density with the constraints derived from the 
observations of cosmic antiprotons and GC 7-rays. This fig- 
ure shows the result for T ro = 10~ 3 and u = 2 (RSII brane 
cosmology scenario §]), where v is defined in Eq. (|3J) for the 
Hubble rate enhancement function. We show the result for 
three different halo models. The labels on the curves show 
the origin of the bound. 



those WIMPs which fulfill it in the standard case. The 
possible enhancement of the Hubble rate can therefore 
be constrained by applying at the same time the relic 
density constraint for the DM and the upper bound on 
the WIMP annihilation cross section as derived from the 
indirect searches for DM. 

The argument above built on the crucial assumption 
that the WIMP annihilation cross section (av) is tem- 
perature independent. If this is not the case then the 
constraint from the indirect searches, which bound the 
present (crv)o, could not be directly combined with the 
constraints on the relic density, which depends on the 
cross section in the early Universe. The relation between 
the cross section in the two epochs should be taken into 
account, a relation which would often lead to a looser 
bound on the enhancement of the Hubble rate. A discus- 
sion of this topic was given in Ref. [l[ , where it was shown 
that modifications are usually not very large, unless some 
specific situations, like e.g. coannihilation, occur. In this 
paper we will show our results only for the case where 
the WIMP annihilation cross section is temperature in- 
dependent. 

Let us derive the bound on the enhancement of the 
Hubble rate from the combination of the constraints on 
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FIG. 6: The same as in Fig. [5] for different values of the 
parameter v: v = 2 (RSII brane cosmology scenario §]), 
v = 1 (kination scenario [5fl), v — and u = — 1 (scalar- 
tensor cosmology scenario [J]). The NFW and Moore DM 
density profiles are shown. 



the WIMP relic density and cross section. To obtain the 
WIMP relic density we solve the Boltzmann equation im- 
plemented with the modified Hubble function, Eq. 
Let us therefore again go through the free parameters of 
the enhancement function. The freeze-out temperature, 
Tf, is determined by the WIMP mass and annihilation 
cross section and is therefore not a free parameter. The 
re-entering temperature, T re , is a free parameter, but 
we showed in Ref. [l| that the bound on the Hubble ex- 
pansion is independent of T rc as long as T rc <C Tf. In 
this paper we set T re = IO -3 GeV, which is always much 
lower than the freeze-out temperature and which is the 
lowest value we can assume not to spoil BBN predic- 
tions. As we have mentioned earlier, the exponent v in 
the enhancement function selects the kind of cosmolog- 
ical model. The only true free parameter is therefore 
77, which is normalized as the enhancement of the Hub- 
ble function at the time where the WIMP freezes out in 
standard cosmology. 

To derive the upper bound on the parameter r\ we use 
for each WIMP mass and halo profile the upper bound 
on the annihilation cross section displayed in Fig. [4] 
For a given cosmological model, determined by the value 
of the v parameter, the upper value of 77 is then found 
where the solution of the Boltzmann equation satisfies 
the upper bound of the density constraint Eq. |T|). The 
upper bound on 77 as a function of the WIMP mass is 
shown in Fig. O for v — 2 and in Fig. [5] also for v = 
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— 1,0, 1. Knowing 77, we can calculate the enhancement 
function A(T) at any time once we choose a cosmological 
model. 

Figs. [5] and H show that 7-ray data may be quite effec- 
tive in constraining the preBBN Hubble rate for heavy 
WIMPs, and nicely complement in this large-mass range 
the antiproton results. In order to set bounds more strin- 
gent than antiprotons, however, a steep density profile 
is required: in the case of a NFW distribution, the 7- 
ray observations of HESS arc able to set limits only for 
WIMPs heavier than 1 TcV. For this mass range, how- 
ever, depending on the actual cosmological evolution, the 
bound can be relevant and much stronger than the an- 
tiproton bound: in the case of kination models, the 7- 
rays predicted for a NFW profile limits the maximal Hub- 
ble rate enhancement to be less than a factor of 5000; for 
scalar-tensor cosmologies the maximal enhancement goes 
down to a factor of 500. On the other hand, for a Moore 
profile the bounds are quite stringent: the maximal en- 
hancement of the Hubble rate in this case is a factor of 
100. In addition, for the Moore profile an enhancement 
is not possible for all the WIMP mass exceeding 1 TcV. 

As an example of how our limits can be further used 
to constrain the basic properties of specific cosmological 
models, let us consider the implications for the brane 
Randall-Sundrum II model [3], which in our notations 
corresponds to v = 2 and 77 = ^/p r (Tf)/(2A). Here p T 
is the radiation energy density and A is the tension of 
the brane, related to the 5-dimensional Planck mass M5 

o Af 6 

by the relation A = j^Tp-- For the Moore profile and 

pi 

WIMP masses m x — O(500 GeV) we see from Fig. [5] 
that 77 < O(10 2 ). This implies A > 2 • 10 2 GeV 4 , which 
corresponds to a lower bound on M5 of M5 > 7 • 10 3 
TcV. This is almost two orders of magnitude better than 
what was found in Ref. [l| by using antiproton data 
for the same value of WIMP mass. It is more stringent 
than what can be obtained from other cosmological tests: 
BBN sets the limit M 5 > 13 TcV (HE^, while Ref. @ 
sets M5 > 600 TcV but for a DM candidate in a specific 
supersymmetric model. Microgravity experiments 
still set the best bound M 5 > 10 5 TeV. 



VI. CONCLUSIONS 

In this paper we have discussed the possibility to use 
7-ray data from the GC to constrain the cosmological 
evolution of the Universe in a phase prior to primordial 
nucleosynthesis, namely around the time of CDM decou- 
pling from the primeval plasma. We extended the argu- 
ments already discussed in a previous paper of ours [l| , 
where instead cosmic-ray antiprotons were used. The ba- 
sic idea is that in a modified cosmological scenario, where 
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the Hubble expansion rate is enhanced with respect to 
the standard case, the CDM decoupling is anticipated 
and therefore the relic abundance of a given DM candi- 
date is enhanced. This implies that the present amount 
of CDM in the Universe may be explained by a WIMP 
which possesses annihilation cross section (much) larger 
than in standard cosmology. This enhanced annihilation 
cross section implies larger fluxes of indirect detection 
signals of CDM, due to the annihilation of relic WIMPs 
in the halo of our Galaxy. 

The stringent bounds on the maximal enhancement of 
the preBBN cosmic evolution, determined by the antipro- 
ton signal from DM annihilation and obtained in Ref. 
[H, have been complemented here by 7-ray searches. We 
have shown that the HESS measurements, which refer to 
relatively large 7-ray energies, are able to set constraints 
for WIMPs heavier than a few hundreds of GeV, depend- 
ing on the actual DM halo profile. These results are 
complementary to those coming from antiprotons, which 
instead are important for WIMPs lighter than a few hun- 
dreds of GeV. In the case of a Moore profile, these bounds 
are very strong, and imply that a WIMP heavier than 
about 1 TeV is not compatible with a cosmological sce- 
nario with enhanced expansion rate prior to BBN. Less 
steep profiles provide less stringent bounds, always for 
heavy WIMPs: the NFW halo bounds the maximal Hub- 
ble rate enhancement to be below a factor between 5 x 10 2 
and 5 x 10 5 , depending of the cosmological scenario. On 
the other hand, an isothermal sphere does not provide 
any relevant limit (better than the antiproton bound) 
for any mass. 7-ray data from the EGRET satellite are 
important for intermediate-mass WIMPs, but only for 
the very steep Moore profile. 

Data from the GLAST satellite-based experiment will 
add relevant information for DM particles in a range of 
masses which goes from 100 GeV to a few hundreds of 
GeV, furthcrly complementing the analysis we have been 
able to perform by using antiprotons (relevant for masses 
below 100-200 GeV) and available 7-rays from the GC 
(relevant for masses above about 500 GeV). 

We can therefore conclude that DM indirect detection 
searches, in addition of being a powerful and important 
tool for studying the DM component of the Universe, may 
also have an important role in constraining the cosmic 
evolution, with an impact on dark energy models, modi- 
fied gravity scenarios and theories of extra-dimensions. 
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